function [f,a]=getFrequencySpectrum(data,dt)
N=length(data);
nfft=2^nextpow2(N);
f=1/dt*(0:(nfft/2))/nfft;
W=fft(data,nfft);
a=abs(W);
a=a(1:(nfft/2+1));
a=a/max(a);
end